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ABSTRACT 

We  consider  the  mathematical  behavior  of  a  viscous.  Incompressible  fluid 
bounded  above  by  an  atmosphere  of  constant  pressure  and  below  by  a  horizontal 
bottom.  After  reviewing  the  existence  and  regularity  theory  for  the  equations 
governing  the  motion,  we  establish  rates  of  decay  for  solutions  near 
equilibrium.  Die  function  describing  the  height  of  the  free  surface  decays 
like  t~^2  ;  the  velocity  field  decays  like  t'1 .  These  estimates  are  shown 
first  for  the  linearization  about  equilibrium  and  then  for  the  full  nonlinear 
problem.  Complete  details  will  be  given  elsewhere. 
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SIGNIFICANCE  AND  EXPLANATION 


One  primary  purpose  of  mathematical  analysis  is  to  describe  the  behavior 
of  solutions  of  equations  representing  physical  phenomena.  The  results  shed 
light  on  the  validity  of  the  mathematical  model  as  a  representation  of  the 
physical  process  and  may  in  turn  offer  new  information  about  the  phenomena. 
For  the  problem  considered  here,  the  motion  near  equilibrium  of  a  viscous 
fluid  bounded  by  a  free  surface,  fundamental  properties  of  existence  and 
regularity  were  established  by  the  first  author  in  earlier  work.  This  paper 
continues  the  study  by  describing  further  the  behavior  of  solutions  over  a 
horizontal  bottom.  Solutions  decay  at  a  fixed  rate  because  of  the  effect  of 
viscosity  in  the  interior  of  the  fluid.  The  slow  rate  of  decay  is  determined 
by  the  behavior  of  long  wavelengths  in  the  linearization  about  equilibrium. 
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LARGE-TIME  BEHAVIOR  OF  VISCOUS  SURFACE  WAVES 
J.  Thomas  Beale*  and  Takaaki  Nishida** 

« 1 .  Introduction . 

Me  are  concerned  with  solutions  global  In  tlae  to  a  free 
surface  problem  of  a  viscous  Incompressible  fluid,  which  Is 
formulated  as  follows:  The  motion  of  the  fluid  Is  governed  by 
the  Navier-Stokes  equations 

u^  +  (u*v)u  -  v^u  +  vp  ■  0 

(1.1)  }  In  fi(t), 

v  •  u  •  0 

where  Q(t)  ■  {x  «  R2,  -b  <  y  <  ij(  t,x)>  Is  the  domain  occupied 
by  the  fluid.  The  free  surface  Sy:  y  -  <7(t,x)  satisfies  the 
kinematic  boundary  condition 

(1.2)  <7t  +  Uji?x  ♦  «2^x  "  u3  "  0  on  Sp* 

The  stress  tensor  satisfies  the  free  boundary  condition: 
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(1.3)  pnj  -  «(»!  ^  + 


1 0*7  -  P*<(1  ♦  l»‘»|2)~1/2*<7)ln1  on  Sy, 


where  n  is  the  outward  normal  to  Sp,  g  la  the  acceleration 


of  gravity  and  p  la  the  nondlaenalonallzed  coefficient  of  sur- 


face  tension.  On  the  bottom  S^:  y  ■  -b  we  have  the  fixed 


boundary  condition 


(1.4) 


u  -  0  on  SB. 


We  assume  here  that  b  la  constant,  although  Theorem  1.1  below 


holds  for  b  *  b(x) . 


We  consider  the  Initial  value  problem  (1.1)— (1.4)  with  the 


data  at  t  ■  0 


(1.5) 


I  n  -  *70(x) 


x  «  R2, 


u  -  uQ(x,y) 


in  a0. 


where  a.  ■  fi(0).  Local  existence  theorems  for  (1.1) -(1.5)  are 


proved  for  both  cases  with  or  without  considering  surface  tension 


([l)f[2]).  The  problem  of  existence  global  In  time  for  (1.1)- 


(1.5)  neglecting  the  surface  tension  {p  m  0)  has  a  difficulty 


which  was  pointed  out  in  [1].  However  If  the  surface  tension  Is 


taken  Into  account,  the  following  global  existence  and  regularity 


theorem  is  proved. 
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PM  1.1  ([2]) 

Let  3  <  r  <  7/2.  Aaauaa  tha  compatibility  co ndltlona  on 


tha  Initial  data: 


’  •  “o  ■  0 


<!•«>  {  “‘Vi.Xj  +  «aOl3,x1>n3>t.n  ’ 


In  n0. 


u0  "  0 


0  on  y  -  q 0(x) . 
on  y  -  *?0(x) . 


Thara  axtata  0Q  >  0  auch  that  It  tha  Initial  data  aatlaf y 


“•7)  "o  '  ’’o'h^R2)  +  <  8° 


than  thara  axtata  a  unlqua  global  aolutlon  q, u,p  of  (1.1)- 
(1.0),  which  aatlaflaa 


(1.8)  q  a  Kr+1/2(R+xR2) ,  u  a  Kr(R+xO(t)),  *p  a  Kr”2(R+xO(t) ) 


further ,  glvan  any  >  0  and  any  k  >  0,  thara  axlata  8^  > 
auch  that  It 


(l.») 


B0  <  *1 


than  tha  aolutlon  bacoaaa  aaooth  tor  t  >  ,  I .a . , 


(1.10)  q  a  Kr+k+1/2( (Tj.-JxR2) ,  u«  Kr+k( (Tlf«*)xO(t) ) 
▼p  «  Kr*k”2( (T, ,«•)  x  0(t)). 


In  particular  the  aolutlon  with  k  a  2  is  classical.  Hsrs 
Hr( • )  is  the  usual  Sobolev  space  with  norm  |  |r  on  the  domain 
Kr((T1,T2)  x  fi( t))  is  cosposed  of  the  restriction  to  the 
fluid  doaain  fl(t)  of  the  functions  belonging  to 

(1.11)  Kr((TlfT2)xR3)  - 

H°((T1,T2);Hr(R3))  n  Hr/2( (Tj .Tj) ,H°(R3) ) . 
q  «  KP(R+  x  R2)  is  defined  as  follows:  <7  •  Kr((0,T)  x  R2)  for 

p  ▲  9 

any  T  >  0  and  n  m  n  ±  +  H  2  *uch  that  7,  «  K  (R  x  R  )  and 
<72  is  the  Fourier  transform  in  space-tine  of  L1  function  of 
bounded  support.  See  [2]  for  the  details  of  the  function  spaces. 

In  this  summary  we  give  an  asymptotic  decay  rate  for  the 
solution  of  the  above  theorem. 

1  2 

Theorem  1.2.  if  rjQ  e  L  (R  ),  then  there  exists  «2  >  0  such 
that  if 

(1.12)  Ej  1  EQ  +  |i70|  x  <  fl2  , 

L 

then  the  aolutlon  has  the  decay  rate: 

l«a*(t)l0  <  CEjd  +  t)”(1+c0/2, 

|u(t)|2,  ( vp(  t )  |  Q  <  CE^l  +  t)"1. 


a  -  0,1,2 


In  t 2  we  transform  the  free  boundary  problem  (1.1) -(1.5)  to  an 
equivalent  one  on  a  fixed  domain  and  reduce  the  components  of  the 
stress  tensor  to  zero.  The  linear  decay  estimate  is  discussed  In 
ft 3  and  the  nonlinear  one  in  ft4. 


« 2 .  Reduction  of  the  Problem. 

We  recall  some  of  the  main  ideas  for  the  reduction  of  the 

free  surface  problem  in  [2].  First  we  use  the  transformation  of 

the  free  boundary  problem  (1.1)-(1.5)  to  that  on  the  fixed 

2 

(equilibrium)  domain:  a  -  (x  e  R  f-b  <  y  <  0} .  Given  q(t,x) 


we  extend  it  for  y  <  0  as  follows: 


(2.1) 


<7 ( t ,x,y)  -  F_1(e,t,y  i(t ,t)). 


where  q(t,t)  is  the  Fourier  transform  with  respect  to  x  and 
F_1  is  the  Inverse.  If  q(t,-)  belongs  to  Hs(Sp),  then 

belongs  to  H  '  (0),  where  S?  now  denotes  the  upper 
surface  y  ■  0  of  a.  For  each  t  >  0  we  define  the  transfor¬ 
mation  9  on  a  onto  a(t)  *  (x  «  R2,-b  <  y  <  *7(t,x)}  by 


-  .V 
w  «.  - 


(2.2) 


• (Xj.Xj.yjt)  «  (x1,x2,q  +  y ( 1  +  q/b)) 


The  vector  field  u  on  8(t)  •  •(«)  is  defined  from  the  vector 
field  v  on  a  by 


•  «  ".V  •  *  V  *  *  *  tt*  .  *  .  '.V  •  . 

“  •  *  •  •  *  *  •/  **  •  *  •  •  • 


%  *  *  ^  »  1 *  <*  •»  •  *  ■  4  •  *  .  *  .  ’ 


Here  we  have  gathered  the  linear  terns  on  the  left  hand  side  and 
all  the  nonlinear  terns  on  the  right  hand  side  of  the  equation. 


Next  we  reduce  the  tangential  component  of  the  atreaa  tensor 

/ O 

rA,  i  *  1.2  to  zero:  Given  F1  •  H  (Sp),  1  -  1,2,  choose 

r+1 

the  vector  z  «  H  (n)  satisfying  the  condition 

2 

j  z  ■  0,  z  ■  0,  z  -  (Fj.-Pj.O)  on  Sp, 

1  z  •  0,  fly  z  ■  0  on  S0. 

Then  w  -  v  x  z  satisfies 


Wi,X3  +  ¥l3,Xi  "  pi' 

i  - 

1,2  on  Sp, 

v  •  w  ■  0 

in 

a  , 

w3  -  0 

on 

V 

Therefore  17,  v'  «  v  -  w,  q  satisfy  the  system  (2. 4) -(2.9)  with 
the  replacements  F  by  F^  ■  F  -  wfc  +  uA*  and  Fi ,  i  »  1 , 2  by 
0.  The  prime  in  v'  is  omitted  hereafter. 

Finally  we  rewrite  the  system  (2. 4) -(2. 9)  with  F  ■  F^, 

F^  *  0,  1  -  1,2,  for  q,v,q  in  the  operator  form.  Let  P  be 
the  projection  on  the  subspace  of  solenoldal  vectors  orthogonal 
to  the  subspace  Gr°  «  (vP:  r  «  H 1(ft),  r  ■  0  on  Sp)  of  H°(n), 
1  •  e . , 

(2.10)  H°  -  PH°  e  Gr°. 

Applying  P  to  (2. 5)  we  have 

(2.11)  vt  ”  pPilv  +  ■  PP4’ 

Here  Pvq  can  be  decomposed  to  three  parts  as  follows: 


where  * 


1  -  1,2,3,  are  defined  by 


(1) 

(2.12) 

We  define 

(2.13) 


Using  these 

(2.14) 

(2.15) 

where  f  (17 , 
with  ■ 


2vv3 (X  ,  «  g 17  -  fi*n  •  ■  p3  on  Sp, 

3 

A  ■  0  In  Q  , 

d  k  ^  ^  *  0  on  S_  . 

Y  B 

A  V  «  -P<dV  +  »K ^ 1 ^  , 


l 


R  v 


9 


•  1 2 ) 

R  ((g  -  pA)q )  -  vk'*'  . 


notations  the  systea  (2.4),  (2.11)  has  the  fora 
i7t  -  R  v  , 

vt  +  Av  +  R*((g  -  pA)q)  ■  f  , 

f,vq)  ■  PF4  -  vn*3* .  We  regard  (2.6),  (2.7),  (2.8) 
1  as  doaaln  conditions  on  A. 


3.  Rates  of  Decay  for  the  Linear  Problem 


We  Investigate  the  decay  rate  of  the  solution  of  the 
inearized  equations 


3.1) 

nt  -  Ru  , 

3.2) 

ut  +  Au  + 

R*( (g  -  pA)n)  * 

0  , 

3.3) 

17(0)  «  I7q, 

u(0)  -  u0  at 

t  *  0. 

These 

are  supplied  with  the 

conditions : 

(3.4) 

v  •  u 

« 

o 

M* 

3 

» 

(3.5) 

u,  v  +  u 

*  »  *3  f 

■  0,  i  -  1,2 

on  Sp, 

(3.6) 


u  *  o  on  Sg. 


Theorem  3.1.  ft  E2  -  |i70l  1  +  l^0 1 5/2  +  lu0lo*  Then  thB  ao/u' 

L 

tion  of  (3.1) -{3. 6)  h as  ttia  decay  rata  (t  fc  1): 


(3.7) 


|oai7(t)|0  S  CqE2  t”(1+a)/2,  0  S  a  S  5/2. 


-1 


I u( t ) | 2  S  CQ  E2  t 

The  theorem  is  proved  in  several  steps. 

Let  *  «  (v  -  (q,u)  :  17  «  H1(Sp),  u  «  PH°(n)},  where 
(P**7)j  *  g(p»< 7)0  +  p(vp,vi 7)Q  is  the  inner  product  of  (S^) 


-9- 


‘  •“»■»*  »**  ,*•  ,*»  '»  “• 


**  .,-4 


kA.— 1 


‘  V  * 


s 


and  sat  W  ■  (v:  17  a  HS^2(SF),  u  •  PH2(fl)  and  u  satisfies 
(3.4),  (3.5),  (3.6)).  Let  us  define  the  operator 


0  R)  [  n 

(3.9)  G  v  *  ,  on  D(G)  « 

-R  (g-fiA)i 7  -AJ  u  , 


and  consider  its  closed  extension  which  will  be  denoted  by  G 
again. 


3.2.  The  operator  G  generates  a  contract  Ion  sami group 


e  on  M,  and  W  c  D(G). 


Consider  the  resolvent  equation: 


(3.10) 


The  resolvent  of  G  can  be  extended  in  to  the  left  half-plane. 


as  shown  in  lemmas  3. 3-3. 5. 


3.3.  For  any  rQ  >  0  there  exists  cQ  >  0  such  that  //* 


A  e  (A  -  o  +  lr,  -cn | t |  <  o  <  r  ,  |r|  >  r  ), 


then  the  solution  of  (3.10)  has  the  estimate 


,  * ,  • « *  <  *  * 


‘  *  ••  .*  •*. 


:: 


i 


_  r,  .  .  f  v*m  imj  y.  ■ 


(3.11) 


|U|2  +  |A||U|0  + 


,A’1rU,5/2 

*  c<m0 


+  l*l5/2  +  1*11*13/2 
*  |h|5/2>  • 


He  treat  the  resolvent  near  A  -  0  In  two  cases  separately:  (1) 
the  supports  of  h(t)»  f(t.y)  belong  to  (|(|  t  tQ);  (11)  the 
supports  belong  to  (|t|  *  t  , } .  Here  *  means  the  Fourier 
transform  with  respect  to  x. 


Lemma  3.4.  For  any  (Q  >  0  there  exists  rQ  >  0  such  that  If 
A  «  { | A |  <  rQ)  end  the  support s  of  h(<),  f(t.Y)  balong  to 
{|t |  <fc  Cq)>  than  tha  rasolvant  aquation  (3.10)  has  the  solution 
(J7,u)  satisfying 


(3.12)  lul 2 ' •* 1 5/2  *  C(lh»5/2  +  |f,0)  * 

Let  6(C)  be  the  Fourier  transform  of  0  with  respect  to  x. 

Leamm  3.5.  mere  exist  Cj  >  0  and  ri'r2'  with  v(*/2b)2  >  r2 
>  r^  >  o,  such  that  If  r^  <  |A|  <  r2 ,  than  (A  -  6(C))  1  exists 
for  |t|  <  (j.  for  small  t  there  is  a  ona-qimanslonal  eigen- 
space  which  Is  analytic  with  raspact  to  t>  Tha  alganvatua  and 
eigenvector  have  the  following  axpanslons. 


-11- 


(3.13) 


a  -  -(gb  /3p) |t |  +  0( | e I  )  » 

17  -  1  +  0(  ( {  |  2 )  , 

{  Uj  -  l(gt j/2v) (y2-b2)  +  0(|t|3),  j  -  1,2, 
u3  -  ( g | e | 2/2p) (y3/3-b2y-2b3/3)  +  0(|e|4). 


vv>, 

W‘ 


By  using  lemmas  3. 2-3. 5  the  decay  estlaate  (3.7)  can  be  proved  by 
the  transformation  of  the  Integral  path  of  the  repreeentatlon 


(3.14)  v(t) 


r°+lr  At  -1 
11m  /  eAt(A-0)  V  dA, 

r-»<*  o-ir  0 


a  >  0 


to  the  left  half-plane. 


14.  Nonlinear  Decay  Estimates. 

The  free  surface  problem  (1.1)- (1.5)  was  reduced  to  the 
following  system  In  12. 


(4.1)  <7t  -  R  u  ■  0, 

(4.2)  ut  +  A  u  +  R*(  (g  -  flA)ij)  m  t  , 

(*•3)  *7(0)  -  iJq,  u(0)  ■  uQ, 

where  f  consists  of  nonlinear  terms  depending  on  7,  u,  vp  and 
their  derivatives. 


>v 


The  initial  value  problem  (4.1) -(4. 3)  haa  a  unique  eolution 
(Theorem  1.1)  which  becomes  smooth  for  t  fc  Tj  >  0.  Namely  we 
know  that  if  BQ  <  then 


(4.4) 


(4.5) 


j  <7(t)  «  C(0,«®;H®/2) ,  i  .e. ,  *{t)  «  C(0,»;H3), 


^  u(t)  «  C(0,*;H2),  and 


i  e  K7(  (!,«•)  x  a).  u(t)  «  K6((l,«)  x  fl). 


r  • 

-V-V.v 

L4U_ 


.-'.-•.-.VJ 


*>  -m 


Let  us  define 


•S.N 


M(i7,u;t) 


(  max  sup  (l+s)(1+a)/2|Dai7(s)  |fl  , 
a-0, 1 , 2  OSsSt 


(4.6) 


sup  (l+s)3/2|D3i(s)|0,  sup  ( 1+s) |u(s) | ,) 
ossst  ossdt 


E(q,u)  -  |I7|-  +  |U| 

K7((l,«.)xfl)  K6 ( ( 1 »*)xfl ) 


v  _v. 


a'.v 
\\V  N*.\’ 


Since  the  linearized  equation  is  solved  by  the  semigroup  e  on 
*  in  %3,  this  solution  satisfies  the  variation  of  constants 
formula : 


(4.7)  v(t)  -  etG  vQ  +  J  e(t~s)Q  f(s)de. 


V*  *\  •*-  W1 


«h.r.  v-<’>  ud  *  -  (?„,u,Tp))- 


•  •Vv ■  •'  •  • 


V  *_*  /  ’j  *,* 

*  •  •  •  »  *•  *  *  I 


/  V  *  V  % V  V-’-%V  N*. 


%  %  • 

>yy-y 


•■'Vv'-”  '''.  -'.**  ••••/  V-'! 

1  s'  ‘  1  ■  '  •  •  >  - 


The  first  tera  of  (4.7) 


alrssdy  os t lsatsd  In  »  3 ,  1 . o . , 


(4.3)  [u°(tj  "  ®tG(u°]  *“■  the  decaY  rata: 

|«ai70(t)i0  d  CqE2  t“(1+a,/2.  0  S  a  S  5/2,  i.e., 

\*>an0lt)\Q  S  C0  E2  t’(l4a)/2,  a  -  0,1,2,  and 
(4.9)  {  |D3i7Q(t)  |Q  S  CQ  E2  t"3/2  , 

|u0(t)|2  *  CQ  E2  t"1  . 


It  Is  oufflclsnt  by  (4.4)  to  prove  tho  decay  rate  for  t  a  2. 
Let  us  decompose  the  second  tera  of  (4.7)  Into  three  parts 


t/2 


(4.10)  j‘{e(t‘*,G[°  ]>ds  -  f  ><*•  ♦  f1  >da  ♦  J  {  )<3 

J0  J0  J  t/2  Jt-1 


(t)  + 


(t)  + 


3J 


(t) 


4.1.  let  E-  •  1 17  |  «  ♦  |u|,  |  Dp  |  n .  We  have  the  eet/- 


sete* 


in0  S  C  ( Eg )  {  (  |  D?  |  2  +  |u|2)(|dJ|2  +  |  u|  2  +  |  Dp  |  Q ) 

♦  \n I  „<|D2u|0  +  ! Dp  I q ) ) • 

L 

(4.11)  |DF|0  S  C(E3)((|dJ|2  +  |u|2)(|Dj|3  +  |u|3  ♦  IDp^) 

♦  1^1  .<|D3U|0  +  |D2p|0)). 

L 

|D2F|0  S  C(Eg )  ( (  I  D<7  I  2  ♦  |u|2)(|D«J|4  ♦  |u|3  «■  IDp^) 

♦  nil  ,  ♦  |D7l  .  ♦  I D2i I , ) ( I D4u| 0  +  |D3p|0)> 


Let  »1#  1  -  1,2,3  on  fl  be  the  extension  of  P^  on  Sy  by 

m 

using  if  for  if .  We  have 

(4.12)  I^Iq  S  C(B^)  (  |  D^r  |  2  +  li  I  „)|u|2 

L 

ot  “  a-1 

D  P. ,  a  m  1,2,3,  have  the  ease  estlaate  as  DP  as  above. 

Therefore  f  has  the  same  estlaate  as  (4.11). 

Lenaa  4.2.  Tor  t  a  2,  we  n av*  for  1  ■  1,2 

(4.13)  M(iy  .  ,U.  ;  t )  S  C,  sup  (1  +  s)3/2|f(s)|n 

1  1  Cttast  u 

SC1C(E3)M(i7,U;t)2. 

Also  we  have  for  1  «  3 

(4.14)  M(i7a,u  ;t)  S  C  sup  (1  +  e)3/2|f(e)|, 

*  t-issst 

1  C2C(E3)M(J7,U;t)(M(9,U;t)  +  E(<7,u))  . 

This  Is  proved  by  Theorea  3.1  and  by  Leaaa  4.1  and  Sobolev- 

Mlrenberg's  Inequality  [4].  In  particular,  we  know  that  since  by 

-  6  5 

(4.5)  ij(t )  and  u(t)  are  bounded  In  B  and  H  respectively 

for  any  t  a  1,  we  have  by  (4.6)  for  any  t  a  1 


|D3+a  q(t)|Q  S  C  t(a"3)/2,  a  -  1,2, 

|D2+“  u( t ) | Q  S  C  t(a/3_1),  a  -  1,2. 

Proof  of  Theorea  1.2.  It  follows  from  (4.9)  and  Leans  4.2  that 


M(«7,u;t)  S  CqE2  +  CjdMq.Ujt)2 

♦  C2CM(i7,u;t)  (M(i7,u;t)  +  E(<7,u)). 

Therefore  there  exists  «2  >  0  such  that  If  E2  <  «2,  EQ  <  0J 
then 

M(i7,u;t)  i  C  E2  . 

This  proves  Theorea  1.2. 


Rsnark .  If  the  fluid  has  an  Infinite  depth,  the  eigenvalue  of 
6(t )  has  the  following  expansion: 


(4.15) 


A(e) 


2v|tr  -  oc 1 1 1 


11/4, 


which  Is  quite  different  froa  (3.13).  The  details  will  be 
published  elsewhere. 


v*. 
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